x=wgn(1,1000,10);
y=sin(2*pi*0.01.*[1:1000]);
plot(1:1000,x+y)